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Abstract 

We examine the soft-gluon resummation effects, including the exact spin correlations among 
the final state particles, in the search of the Standard Model Higgs boson, via the process 
gg ^ H ^ WW/ZZ — > 41eptons, at the Tevatron and the LHC. A comparison between the 
resummation and the Next-to-Leading order (NLO) calculation is performed after imposing various 
kinematics cuts suggested in the literature for the Higgs boson search. For the H ZZ mode, the 
resummation effects increase the acceptance of the signal events by about 25%, as compared to the 
NLO prediction, and dramatically alter various kinematics distributions of the final state leptons. 
For the H WW mode, the acceptance rates of the signal events predicted by the resummation 
and NLO calculations are almost the same, but some of the predicted kinematical distributions are 
quite different. Thus, to precisely determine the properties of the Higgs boson at hadron colliders, 
the soft-gluon resummation effects have to be taken into account. 
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I. INTRODUCTION 



Although Standard Model (SM) explains successfully all current high energy physics 
experimental data, the mechanism of electoweak spontaneous symmetry breaking, arising 
from the Higgs mechanism, has not yet been tested directly. Therefore, searching for the 
Higgs boson (H) is one of the most important tasks at the current and future high energy 
physics experiments. The negative result of direct search at the LEP2, via the Higgsstrahlung 
process e~^e~ — > ZH, poses a lower bound of 114.1 GeV on the SM Higgs boson mass 
(Mh) On the other hand, global fits to electroweak observables prefer Mh ^ 200 GeV at 
the 95% confidence level j^, while the triviality arguments put an upper bound ~ 1 TeV j^. 

There is currently an active experimental program at the Tevatron to directly search for 
the Higgs boson. The Large Hadron Collider (LHC) at CERN, scheduled to operate in late 
2007, is expected to establish the existence of Higgs boson if the SM is truly realized in 
Nature. At the LHC, the SM Higgs boson is mainly produced through gluon-gluon fusion 
process induced by a heavy (top) quark loop. Once being produced, it will decay into a 
fermion pair or vector boson pair. The strategy of searching for the Higgs boson depends 
on how it decays and how large the decay branching ratio is. If the Higgs boson is lighter 
than 130 GeV, it mainly decays into a bottom quark pair (bb). Unfortunately, it is very 
difficult to search for the Higgs boson in this mode due to the extremely large Quantum 
Chromodynamics (QCD) background at the LHC. However, the if — > 77 mode can be used 
to detect a Higgs boson with the mass below 150 GeV [J, 15| though the decay branching 
ratio of this mode is quite small, ~ O(10~'^). If the Higgs boson mass (Mh) is in the region 
of 130 GeV to 2Mz {Mz being the mass of Z boson), the H ZZ* mode is very useful 
because of its clean collider signature of four isolated charged leptons. The H WW^*^ 
mode is also important in this mass region because of its large decay branching ratio. When 
Mh > 2Mz, the decay mode H ZZ i~^£~i'^£'~ is considered as the "gold-plated" mode 
which is the most reliable way to detect the Higgs boson up to Mh ~ 600 GeV because the 
backgrounds are known rather precisely and the two on-shell Z bosons could be reconstructed 
experimentally. For Mh > 600 GeV, one can detect the H ZZ £^£~i'u decay channel 
in which the signal appears as a Jacobian peak in the missing transverse energy spectrum. 

The discovery of the Higgs boson relies on how well we understand the signals and its 
backgrounds, because one needs to impose optimal kinematics cuts to suppress the huge 
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backgrounds and enhance the signal to background ratio {S/B). Many works have been 
done in the hterature to calculate the higher^order QCD corrections to the (dominant pro- 

H,y. In 



duction process of the Higgs boson gg ^ H 



BflflflQInl, [i2l,Q,Q: 
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addition to determine the inclusive production rate of the Higgs boson, an accurate predic- 
tion of kinematics of the Higgs boson is very essential for the Higgs boson search. However, 
a fix order calculation cannot reliably predict the transverse momentum (Qt) distribution 
of Higgs boson for the low Qt region where the bulk events accumulate. This is because 
of large corrections of the form \ia{Q^/Q'^) due to non-complete cancellations of soft and 
collinear singularities between virtual and real contributions, where Q is the invariant mass 
of the Higgs boson. Therefore, one needs to take into account the effects of the initial state 
multiple soft-gluon emissions in order to make a reliable prediction on the kinematic distri- 



butions of the Higgs boson. One approach to achieve this is to include parton showering 



which resums t 
and PYTHIA 



le universal leading logs in Monte Carlo event generators, e.g. HERWIG [20 



19| 



2ll |. which are commonly used by experimentalists. The showering process 



just depends on the initial state parton and the scale of the hard process being considered. 
The advantage is that it could be incorporated into various physics processes. Recently, 



an approach to match N 
generators, MC@NLO [2^ 



lO matrix element calculation and parton showing Monte Carlo 

n 

23|, has been proposed. Another approach is to include correctly 



the soft-gluon effects is to calculate an analytical result by using the Collins-Soper-Sterman 



(CSS) resummation formalism 



24, 
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26l . l27| to resum these large logarithmic corrections to 



all order in as- However, in practice the power of logarithms included in Sudakov exponent 



depends on which level the fixed order calculation has been performed 
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32|. 



It is very interesting to compare the predictions between parton showe ring and resumma- 
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tion calculation and detailed comparisons have been presented in Ref. 
which concluded that all of the distributions are basically consistent with each other, except 
PYTHIA in the small Qt region and HERWIG in the large Qt region. 

In addition, the spin correlation among the Higgs decay products has been proved to 



be crucial to suppress the backgrounds 
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381 . Hence, an accurate theoretical prediction. 



which incorporates the initial state soft-gluon resummation effects and the spin correlations 
among the Higgs decay products, is needed. In this paper, we present such a calculation 
and study the soft-gluon resummation (RES) effects on various kinematics distributions 
of final state particles. Furthermore, we examine the impact of the RES effects on the 
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Figure 1: Tree level Feynman diagram of process gg ^ H ^ > ^i^2)V2(— > ^3^4). 



acceptance rate o 



the literature 



371, 



' the signal events with various kinematics cuts (which were suggested in 
, 3^ for Higgs search) and compare them with the leading order (LO) and 
NLO predictions ^. 

The paper is organized as follows. In Sec. [Tll we present our analytical formalism of 
the CSS resummation. In Sec. IIIIl we present the inclusive cross section of the signal 
process for several benchmark masses of the Higgs boson. In Sec. HVl we study the process 
gg^ H ^ WW^*'^ £+£'~ueUi' for Mh = 140 GeV at the Fermilab Tevatron and for Mh = 
170 GeV at the LHC. In Sec. M, we examine the process gg ^ H ^ ZZ^*^ i+ £-£'+£'- for 
Mh = 140 GeV and 200 GeV, respectively, and the process gg ^ H ^ ZZ i^i^uu for 



Mh = 600 GeV, at the LHC. Our conclusions are given in Sec. 



II. TRANSVERSE MOMENTUM RESUMMATION FORMALISM 

At the hadron colliders, the SM Higgs boson is mainly produced via gluon-gluon fusion 
process through a heavy quark triangle loop diagram, cf. Fig [H, in which the effect of the 
triangle loop is replaced by the effective ggH coupling (denoted as the bold dot). Taking 
advantage of the narrow width of the Higgs boson, we can factorize the Higgs boson pro- 
duction from its sequential decay. The resummation formula was already presented in Ref. 



^ The NLO Quantum Electrodynamics (QED) and electroweak (EW) corrections to the Higgs decay process 
H WW/ZZ ^ Ai were calculated in Ref. and Ref. j4l|, respectively Recently, the NLO QCD 



correction to the Higg s boson decays H — *■ WW/ZZ 4g with hadronic four-fermion final states was 
calculated in Ref. [42|- Since the higher order corrections for Higgs production are dominated by the 



initial state soft-gluon resummation effects, we focus our attention on the RES effects in this work. It is 
worth mentioning that the NLO QED corrections to the Higgs boson decay H — ^ WW/ZZ 41 have 
been implemented in ResBos j^sl program, and the phenomenological study of the combined RES effects 
and the QED correction will be presented elsewhere. 



4 



28|. Here, we list some of the relevant formulas as follows, for completeness: 

dQ'^dQ'^dyd(f)Hd^A 



ao{gg H) 



X 



X 



M{H ^ V1V2 ^ ^l^2hU) 



(27r 



where Q, Qt, y, and 4>h are the invariant mass, transverse momentum, rapidity, and az- 
imuthal angle of the Higgs boson, respectively, defined in the lab frame, and dU^ represents 
the four-body phase space of the Higgs boson decay, defined in the CoUin-Soper frame j^^. 
In Eq. ([T]), ■ ■ )|^ denotes the matrix element square of the Higgs boson decay and reads 
as 



MiH^ V1V2 ^ W^a) 



where ray is the vector boson mass, qi{pi) denotes the momentum^ of the vector boson Vi 
(the lepton and Gp is the Fermi coupling constant. Here, 

C± = (^12 + ^12) («34 + ^34) ± 4ai2&12a34&34, 

where 012 and 612 respectively denote the vector and axial vector components of the ^£1^2 
coupling, while 034 and 634 are the ones for Vi^i^. For the W boson, my = mw, and 

a = b=V2, 

while for the Z boson, my = niz, and 

a =4sine2^-l, b = -1 for Z^ti', 
a = 1, b = 1 for Z — > uu, 



The direction of momentum pi is defined to be outgoing from the mother particle. 
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where 6w is the weak mixing angle. In Eq. ([T]), the function Wgg sums over the soft gluon 
contributions that grow as x [lor hi^Q'ip/Q'^)] to all order in ag, which contains the 
singular part as Qt — > 0. The contribution which is less singular than those included in Wgg 
is calculated order-by-order in as and is included in the Y term. Therefore, we can obtain 
the NLO results by expending the above resummation formula, i.e. Eq. ([T]), to the a| order. 



More details can be found in Ref. 



In our calculation, ctq includes the complete LO 



contribution with finite quark mass effects 
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48l |. It has been shown that this 



prescription approximates well the exact NLO inclusive Higgs production rate. 

For the numerical evaluation, we chose the following set of SM input parameters [49||: 

Gf = 1.16637 X 10~^GeV-^ a = 1/137.0359895, 
mz = 91.1875 GeV, a,(m^) = 0.1186, 



rrif, 



0.5109997MeV, = 0.105658389 GeV. 



Following Ref. [5^, we derive the W boson mass as mw = 80.385 GeV. Thus, the square 
of the weak gauge coupling is g'^ = A\/2'm'^GF- Including the 0{as) QCD corrections to 
W qq', we obtain the W boson width as Tw = 2.093 GeV and the decay branching 
ratio of Br(Vr — > iu) = 0.108 5JJ- In order to include the effects of the higher order 
electroweak corrections, we also adapt the effective Born approximation in the calculation of 
the H ZZ — > 4 leptons mode by replacing the sin^ 9\y in the Zii coupling by the effective 
sin^6'^''^ = 0.2314, calculated at the niz scale. 



III. INCLUSIVE CROSS SECTIONS 



For the mass of the Higgs boson being within the intermediate mass range, it will princi- 
pally decay into two vector bosons which sequentially decay into either lepton or quark pairs. 
Leptons are the objects which can be easily identified in the final state, so the di-lepton decay 
mode is regarded as the "golden channel" due to its clean signature and well-known back- 
ground. The drawback is that the di-lepton mode suffers from the small decay branching 
ratio for the vector boson decay {V — > il). For example, the branching ratio of Z ^ £~^i~ 
is only about 3.4%. Due to the huge QCD backgrounds, the purely hadronic decay modes 
are not as useful for detecting the Higgs boson. 

In this paper, we focus on the purely leptonic decays of the vector bosons in the H —>■ 
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Table I: Inclusive cross sections oi gg ^ H ^ VV 41 at the Tevatron Run 2 and the LHC in 
the unit of fb, i.e. a{gg H) x Br{H VV) x Br{V ^1^2) x Br{V isQ for various Higgs 
boson masses. Here, i and i' denote either e or ji. 
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lyVI^^*^ and H ZZ^*^ modes. To cover the intermediate mass range, we consider the 
following benchmark cases: (i) H WW^*^ t^t'-VfVf,, {£, f = e or /i) for Mh = 140 GeV 
at the Femilab Tevatron Run 2 (a 1.96 TeV pp collider), and for Mh = 170 GeV at the LHC 
(a 14 TeV pp collider); (ii) H ZZ^*^ i+i-£'+i'~ {£,£' = e or /i) for Mh = 140 and 
200 GeV at the LHC; (iii) H ^ ZZ ^ t^tvv for Mh = 600 GeV at the LHC, where £ = e 
or /i, and = z/g, z/^ or z/,-. All the numerical results are calculated by using ResBos 

H- We 

adapt CTEQ6.1L parton distribution function in the LO calculation and CTEQ6.1M parton 



distribution function [52| in the NLO and RES calculations. The renormalization scale (fin) 
and factorization scale (/i^) are chosen to be the Higgs boson mass in our calculations, i.e. 
Aii? = /iF = Mh- 

The inclusive cross sections for those benchmark masses of the Higgs boson are sum- 
marized in Table [T] where different searching channels are considered. For comparison, we 
show the Qt distributions calculated by using the RES and NLO calculations in Fig. [3](a). 
The RES calculation is similar to that presented in Ref. [28|, 
57|, but with A^^ included, where 

11C(3) 



B 
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29| with the known A and 
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where Ca = 3, Cp = 4/3, Nf = 5 and the Riemann constant C(3) = 1.202... . We also 
use the modified parton momentum fractions xi and X2 to take into account the kinematic 



corrections due to the emitted soft gluons 



28|, with Xi = mr^^ l^~S and X2 = rriTe ^/VS, 
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Figure 2: Normalized distributions of transverse momentum of Higgs boson predicted by RES 
calculation at the LHC. 



where = \/ Q\ + Q"^ and is the center-of-mass energy of the hadron collider. We 
also adopt the matching procedure described in the Ref. |43| and the non-perturbation con- 
tribution ly^^ of BLNY form in the Ref. Ql- 

In Fig. [21 we show the transverse momentum 
distributions of Higgs boson predicted by RES calculation at the LHC. As we see that the 
peak position is shifted to larger Qt region and the shape becomes broader when the mass 
of Higgs becomes heavier. 

It is clear that the prediction of NLO calculation blows up in the Qt region and the 
RES effects have to be included to make a reliable prediction on event shape distributions. 
In the NLO calculation, it is ambiguous to treat the singularity of the Qt distribution near 
Qt = 0, see the dashed curve in Fig. [3](a). Before presenting our numerical results, we 
shall explain how we deal with the singularity in the NLO calculation when Qt ~ 0. In 
ResBos, we divide the Qt phase space with a separation scale Qt^. We calculate the Qt 
singular part of real emission and virtual correction diagrams analytically and integrate the 
sum of these two parts up to Q'^r^^ . By this procedure, it yields a finite NLO cross section, 
for integrating Qt from up to Q^^, which is put into the Qt = bin of the NLO Qt 
distribution (for bin width larger than Qt^). Since the separation scale Q^^^ is introduced 
in the theoretical calculation for technical reasons only and is not a physical observable, the 
sum of both contributions from Qt > Qt^ and Qt < Qt^ should not depend on Qt^. As 
shown in Fig. [3](b), the NLO total cross section indeed does not depend on the choice of 
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Figure 3: (a) Distribution of transverse momentum of Higgs boson, and (b) NLO total production 
cross section of Higgs boson via gluon gluon fusion as Mu = 170 GeV at the LHC. 

Q^-*^^ as long as it is not too large. We refer the readers to the Sec. 3 and the Appendix 



of Ref. |43!| for more details. In this study, we choose Q^-^^ = 0.96 GeV in our numerical 
calculations. 

As mentioned in the Introduction, MC@NLO, which matches NLO calculations and par- 
ton showering Monte Carlo event generators, not only predicts a reliable Qt of the Higgs 
boson but also includes spin correlations among the Higgs decay products. Therefore it is in- 
teresting to compare the Qt predictions between MC@NLO and RES calculations. In order 
to compare the differences in shape more precisely, we show the Qt distributions predicted 
by MC@NLO and ResBos in Fig. H for Mh = 140(170, 200, 600) GeV. All distributions 
are normalized by the total cross sections for the corresponding Higgs boson masses. The 
bottom part of each Qt distribution plot presents the ratio between MC@NLO and ResBos. 
We note that for a light Higgs boson the distributions are consistent in the peak region 
34 l36|, where the difference is about 10%, but they are quite different in the large Qt 
region, say Qt > 100 GeV. For a heavy Higgs boson, e.g. Mh = 600 GeV, these two dis- 
tributions are very different in the small Qt region, and MC@NLO tends to populate more 
events in the small Qt region, as compared to ResBos. Since the Higgs boson is a scalar, the 
distributions of Higgs boson decay products just depend upon the Higgs boson's kinematics. 
Therefore, the difference in the Qt distribution predictions between MC@NLO and ResBos 
may prove to be crucial for the precision measurements of the Higgs boson's properties. A 
further detailed study of the impact of the Qt difference on the Higgs boson search is in 



9 




1 — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
= 140 GeV (a) 



- RES 
-- MC@NLO 

-~1 -1 I I ' ' ' ' L. 



"T — I — I — I — r 



1.5 
1.0 
0.5 

0.03 
0.02 
0.01 



1.5 
1.0 
0.5 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 
LMC@NLO / RES,^^'-~\ 


1 1 1 1 1 1 - 


- 


- 



50 100 150 200 

Qt (GeV) 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



Mj^ = 200 GeV (c) 




J I I I III I I 1 1' I I I I I 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I I I I I I I I L 



50 100 150 200 

Qt (GeV) 



1 — I — I — I — I — I — I — I — I — r 
Mj^=170GeV (b) 




50 100 150 200 



Qt (GeV) 




50 100 150 200 

Qt (GeV) 



Figure 4: Comparison of the Qt distributions between ResBos and MC@NLO. 
order and will be presented elsewhere. 



IV. PHENOMENOLOGICAL STUDY OF THE H WW MODE 

In the search for SM-like Higgs boson via H WW^*^ mode, two scenarios of W boson 



decay were considered in the literature 



60 



61 



one is that both W bosons decay 



leptonically, another is that one W boson decays leptonically and another W boson decays 
hadronically. Throughout this paper, we only concentrate on the di-lepton decay mode, i.e. 
H — > l/^VT*^*'' — > i'~ uiUii , at the Tevatron and the LHC. The collider signature, therefore, 
is two isolated opposite-sign charged leptons plus large missing transverse energy {^t) which 
originates from the two neutrinos. In this section, we first examine the RES effects on various 
kinematics distributions, and then show the RES effects on the Higgs mass measurement. 
Finally, we study the RES effects on the acceptances of the kinematics cuts suggested in the 
literature for Higgs search. 
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Figure 5: Kinematic configurations of Higgs decay [H WW — > i~^£ vv) in the rest frame of H: 
(a) H W^W^, (b) H WlWZ and (c) H W^W^ . Here, +(-, 0) denotes the right-handed 
(left-handed, longitudinal) polarization state of the W boson. The long arrows denote the moving 
directions of the final-state leptons. The short bold arrows denote the particles' spin directions. 

A. Basic kinematics distributions 



For a heavy Higgs boson, the two vector bosons, which are generated from the spin-0 Higgs 
boson decay, are predominantly longitudinally polarized, while the longitudinal and trans- 
verse polarization states are democratically populated when the Higgs boson mass is near the 
threshold for decaying into the vector boson pair [63|, l6j|. When 140 GeV < < 170 GeV, 
the transverse polarization modes contribute largely. The two charged leptons in the final 
state have different kinematics because of the conservation of angular momentum, cf. Fig. 
[HI therefore, one charged lepton is largely boosted and its momentum becomes harder while 
another becomes softer. Making use of these differences, one can impose asymmetric trans- 
verse momentum (p^) cuts on the two charged leptons to suppress the background. On the 
event-by-event basis, we arrange the two charged leptons in the order of transverse momen- 
tum: p^"""" denotes the larger pr between the two charged leptons while is the smaller 
one. Fig. [6] shows the distributions of p^"""", -p^ and missing energy for Mh = 140 GeV 
at the Tevatron (first row) and for Mh = 170 GeV at the LHC (second row). Furthermore, 
in Fig. [7] we show the distributions of cos9ll, (pn and AY/^^. without imposing any kinemat- 
ics cut, where cosOll is the cosine of the opening angle between the two charged leptons, 
(f)LL is the azimuthal angle difference between the two charged leptons on the transverse 
plane, and AYll is the rapidity diff'erence of two charged leptons in the lab frame. Since 
we are mainly interested in the shapes of the kinematics distributions, the curves shown in 
the figures are all normalized by the corresponding total cross sections. The solid curves 
present the distributions including the RES effects, the dashed and dotted curves present 
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Figure 6: Normalized distributions of the leading transverse momentum prp'^""^ , softer transverse 
momentum of the leptons, and the missing energy 99 ^ H ^ WW — > (.'^ u^vii . The 

panels (a) to (c) are for Mh = 140 GeV at the Tevatron, and (d) to (f) are for Mh = 170 GeV at 
the LHC. 

the distributions calculated at the NLO and LO, respectively. 

We note that the pt distributions of the charged leptons and the missing energy distri- 
butions are modified largely by the RES effects. This can be understood as follows. The 
two charged leptons prefer to move in the same direction due to the spin correlation among 
the decay products of the Higgs boson, cf. the distributions of cos 6'/,^ in Figs.[7](a) and (d). 
Hence, one can approximately treat the Higgs boson decay as "two-body" decay, i.e. decay- 
ing into two clusters as H (z/^z/^/). This is in analogy to the W boson production 
and decay in the Drell-Yan process, ud — » d^u^ which has been shown in Ref. [51 1 
that the transverse momentum of lepton (p^) is very sensitive to the transverse momentum 
of the W boson. The same sensitivity also applies to I^t- As shown in Figs.[6](c) and (f), the 
clear Jacobian peak of the distribution around Mh/'2. in the LO calculation is smeared 
in the NLO and RES calculations. Furthermore, the distribution in the NLO and RES 
calculations has a long tail due to the non-zero transverse momentum of the Higgs boson. 
Since the RES calculation includes the effects from multiple soft-gluon radiation, the 
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Figure 7: Normalized distributions of cos Oil, 4>ll and ^Y^l va. gg ^ H ^ WW — > vivur. 
The panels (a) to (c) are for Mu = 140 GeV at the Tevatron and (d) to (f ) are for Mh = 170 GeV 
at the LHC. 

distribution near the Jacobian peak is further smeared in the RES calculation as compared 
to the NLO calculation. When Mh = 140 GeV, only one W boson is on-shell and the two 
charged leptons do not move as close as they do in the case of Mh = 170 GeV (in which 
case, both W bosons are on-shell). However the parallel configuration is still preferred. 

The dominant backgrounds of the H WW^*^ mode are from the W boson pair pro- 
duction and top quark pair production. The latter, as the reducible background, can be 
suppressed with suitable cuts such as jet-veto, but the former, as the irreducible back- 
ground, still remains even after imposing the basic kinematic cuts. In order to reduce this 
intrinsic background, one needs to take advantage of the characteristic spin correlations of 
the charged leptons in the H WW^*'^ — > £~^£'~U£i'ii decay. For example, the distribution of 
the difference in azimuthal angles of the charged leptons peaks at smaller value (cf. Figs.lWb) 
and (e)) for the signal than that for the WW continuum production background [60|, 
We note that the RES effects do not affect the cos 9ll and (f>LL distributions very much, as 
shown in Figs.[7](a), (b), (d) and (e). 

To closely examine the difference in their predictions, we also present the ratio of the RES 
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Figure 8: Ratio of the Resummation contribution to NLO and LO contributions m. gg ^ H —>■ 
WW^*^ l-^i'-ueui'. The panels (a) to (c) are for Mh = 140 GeV at the Tevatron while (d) to (f) 
are for Mh = 170 GeV at the LHC. 

contribution to the NLO and LO contributions in Fig.[8l We note that the ratio is about one 
below the peak regions of py™""", and ^t, and becomes larger than one above the peak 
region, where both the LO and NLO contributions drop faster than the RES contribution 
does, which is consistent with the results shown in Fig. [6l This uneven behavior indicates 
that one cannot simply use the leading order kinematics with the constant i^'-factor included 
to mimic the higher order quantum corrections. We should stress that even though the NLO 
and RES calculations include the same contributions of the hard gluon radiation from initial 
states, the effects of the multiple soft-gluon radiation could cause more than 25% difference 
between RES and NLO predictions in the large px and region. 



B. Higgs mass measurement 

In order to identify the signal events clearly, it is crucial to reconstruct the invariant mass 
of the Higgs boson. Unfortunately, one cannot directly reconstruct the Mh distribution in 
the H — i> WW mode due to the two neutrinos in the final state. Instead, both the transverse 
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mass Mt and the cluster transverse mass M( 



c 



65|, defined as 



Mt = j2p^^^T(l-cosA0(pf^^,^T)), 



Mc = Vp^^' + ml^+ ^T, (3) 

yield a broad peak near Mh- In Eq. (JSj), pff^ {rriLL) denotes the transverse momentum 
(invariant mass) of the two charged lepton system, and A(f){pff^, ^t) is the diff'erence in 
azimuthal angles between and on the transverse plane. We note that the upper 
endpoint of Mt distribution can clearly reflect the mass of Higgs boson, cf. Figs. [9](a) and 
(c). Mt is insensitive to Qt because it depends on Qt in the second order, cf. Eq. ([3]). 
Therefore, the position of the endpoint is only subject to Mh and Vh- The latter effects can 
be safely ignored because is very small (less than about 1.5 GeV), for the Higgs boson 
mass less than 200 GeV. The cluster transverse mass Mc also exhibits a clear Jacobian peak 
with a clear edge at Mhj cf. Figs. [9](b) and (d). But both the line shape and the Jacobian 
peak of Mc distribution are modifled by the RES efl'ects because Mc is directly related to 
I^T which depends on Qt in the first order. We suggest that one should use Mt to extract 
the mass of Higgs in H WW^*^ mode because the upper endpoint of the 

Mt distribution is insensitive to high order corrections. 



C. Acceptance study 

In order to separate the signal from its copious backgrounds, one needs to impose optimal 
cuts to suppress backgrounds and enhance the signal to background ratio {S/B ) simulta- 
neously. The selection of the optimal cuts highly depends on how well we understand the 
kinematics of the signal and background processes. As shown above, the RES effects modify 
the distributions of transverse momentum of the charged leptons and the missing energy 
largely, therefore, it is important to study the RES effects on the acceptances of the kine- 
matics cuts. Here, we impose a set of kinematics cuts used by experimental colleagues in 



Refs. 



37 



39||. The corresponding acceptances are summarized in Table | 



For the search for a 140 GeV Higgs boson at the Tevatron, we impose the following 
basic cuts: 

pLma. ^ 15 GeV , > 10 GeV, 

\Yl\ < 2.0 , ^T > 20 GeV, (4) 
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Figure 9: Normalized distributions of the transverse mass Mt and the cluster mass Mc in gg — > 
H WW^*'^ l+i'-VfDe,! : (a) and (b) are for Mh = 140 GeV at the Fermilab Tevatron while (c) 
and (d) are for Mh = 170 GeV at the LHC. 

and the optimal cuts as follows: 



where Yi denotes the rapidity of charged lepton, and Ht denotes the scalar sum of 
the transverse momenta of final state particles, i.e. Ht = \Pt\ + \Pt\ ~^ \ Pt\- The 
overall efficiency of the cuts is about 68% , 69% and 70% after imposing the basic cuts 
(Eq. ([11)) for RES, NLO and LO calculations, respectively, and about 44% for both 
RES and NLO calculations and 46% for LO calculation after imposing the optimal 
cuts (Eq. 

• For the search of a 170 GeV Higgs boson at the LHC, we require the following basic 
cuts : 




+ 20 GeV <HT<Mi 



< Mt < Mh - 10 GeV 



(5) 



pL^a. > 20 GeV , pff> 10 GeV, 



Yl\ < 2.5 , ^T > 40 GeV, 



(6) 
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Table II: Acceptance oi gg ^ H ^ WW'^*^ — > I'^l' uiv^i events after imposing the basic cuts and 
the optimal cuts for Mh = 140 GeV at the Tevatron and Mh = 170 GeV at the LHC. 





Mh = 


140 GeV 


Mh = 


170 GeV 




basic (Eq. (gj) 


optimal (Eq. ^) 


basic (Eq. ^) 


optimal (Eq. 0) 


RES 


0.68 


0.44 


0.61 


0.19 


NLO 


0.69 


0.44 


0.61 


0.19 


LO 


0.70 


0.46 


0.63 


0.20 



and the optimal cuts: 

rriLL < 80.0 GeV , Mh - 30.0 GeV < Mt < Mh , 

(j)LL < 1.0 rad , Oll < 0.9 rad , \AYll\ < 1.5 , (7) 

The cut efficiency is about 61% for both RES and NLO calculations, but about 63% 
for LO contribution after imposing the basic cut (Eq. IQ). After imposing the optimal 
cuts (Eq. ([71)), the acceptances of RES and NLO are about 19%, while LO is 20%. 



V. PHENOMENOLOGICAL STUDY OF THE H ^ ZZ MODE 

In the search for the SM Higgs boson, the H ZZ^*^ is an important 

discovery channel for a wide range of Higgs boson mass. The appearance of four charged 
leptons with large transverse momenta is an attractive experimental signature. This so- 
called "gold-plated" mode provides not only a clean signature to verify the existence of the 



Higgs boson but also an excellent process to explore its spin and CP properties [66||. In 
this section, we study three mass values of Mh (140, 200 and 600 GeV) at the LHC. For 
Mh = 140 GeV and 200 GeV, we require the two Z bosons both decay into charged leptons; 
for Mh = 600 GeV, we require one Z boson decays into a charged lepton pair and another Z 
boson decays into a neutrino pair, i.e. £~^i~i'u. In this section we first study the RES effects 
on various kinematics distributions and then examine the RES effects on the acceptances of 
the kinematics cuts. 
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Figure 10: Normalized distributions of p^™"" and in gg ^ H ^ ZZ ^ e+r£'+£'~: (a) and (b) 
are for Mh = 140 GeV; (c) and (d) are for Mh = 200 GeV at the LHC. 

A. gg^H ^ ZZ^*^ £+i~£'+i'- 



H ZZ^*^ mode in the order of transverse momentum. We denote p;^™""' as the largest 
Pt of the four charged leptons while p-^ the second leading pt- In Fig. [TOl we show the 
distributions of p^""""" and for Mh = 140 and 200 GeV, respectively. Due to the similar 
kinematics discussed in the H WW^*'^ mode, the shapes of the distributions of p^™""' 
and are changed significantly by the RES effects. The typical feature is that the RES 
effects shift the pr of the charged lepton to the larger px region and, therefore, increase the 
acceptances of the kinematics cuts. The numerical results will be shown later. 

Although one can measure the Higgs boson mass by reconstructing the invariant mass of 
the four charged leptons, one still needs to reconstruct the Z bosons in order to suppress 
the backgrounds. The reconstruction of the Z boson depends on the lepton fiavors in the 
final state. In this study, we consider two scenarios: different fiavor charged lepton pairs, 
i.e. H 2e2yU, and four same fiavor charged leptons, i.e. H 4e(or 4/i). Hence, we have 
two methods for reconstructing the Z bosons: 



Similar to the H — > mode, we also arrange the four charged leptons of the 
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1. Different flavor charged lepton pairs (2e2/x): 

In this case, it is easy to reconstruct the Z bosons because both electron and muon 
lepton flavors can be tagged. Using the flavor information, the Z bosons can be 
reconstructed by summing over the same flavor opposite-sign leptons in the final state. 

2. Four same flavor charged leptons (4e/4/i): 

If the flavors of four leptons are all the same, one needs to pursue some algorithms 
to reconstruct the Z boson mass. In our analysis, we flrst pair up the leptons with 
opposite charge. We require the pair whose invariant mass is closest to Mz to be 
the one generated from the on-shell Z boson, and the other pair is the one generated 
from another Z boson, which could be on-sell or ofl'-shell. We name it as the minimal 
deviation algorithm (MDA) in this paper. 

In Fig. dH we show the distributions of the reconstructed Z boson for 140 and 200 GeV, 
respectively. When the flnal state lepton flavors are different, one can reconstructed the Z 
boson perfectly by matching the lepton flavor. For the same flavor leptons, the reconstructed 
Z boson distributions in the MDA are shown as the solid, dashed and dot-dashed curves for 
RES, NLO and LO, respectively. Some points are worthy to point out as follow: 

• We note that the MDA can perfectly reconstruct the distributions of true Z bosons, 
irregardless whether these two Z bosons are both on-shell or only one of them is 
on-shell. 

• When Mu = 200 GeV, both Z bosons are produced on-shell and boosted. The peak 



position of the transverse of momentum is around y (Mh/^)"^ — Tn^z ~ ^1 GeV. For 
all the cases, the RES effects change the shape of largely and shift the to the 
larger value region. 



It has been shown in Ref. [67| that angular correlation between the two Z bosons from 
the Higgs decay can be used to suppress the intrinsic background from ZZ pair production 
efficiently. One of the useful angular variables is the polar angle {6*z) of the (back-to-back) 



Z boson momenta in the rest frame of the Higgs boson 67|. As shown in Fig. [I2l in the rest 
frame of Higgs boson, the back-to-back Z bosons like to lie in the direction perpendicular to 
the 2;— axis, which is the moving direction of the Higgs boson in the lab frame. After being 
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Figure 11: Normalized transverse momentum of Z boson in gg H ZZ^*^ £+£-£'+£'- at the 
LHC: (a) is the px distributions of off-shell Z boson for Mh = 140 GeV, (b) is the pr distributions 
of on-shell Z boson for Mh = 140 GeV and (c) is the pr distributions of on-shell Z boson for 
Mh = 200 GeV. 



T~l — I I I I I I I — I I I I I — ^ 




Figure 12: Normalized polar angle of the (back-to-back) Z boson momenta distributions in the rest 
frame of the Higgs boson in gg ^ H ^ ZZ^*'^ £+£-£'+£'- at the LHC: (a) is for Mh = 140 GeV 
, (b) is for Mh = 200 GeV. 

boosted to the lab frame, two Z bosons will move close to each other, c.f. Fig. [T3l(b). where 
9zz is the opening angle between the two Z bosons in the lab frame. Another interesting 
angular variable is the angle between the two on-shell Z boson decay planes {(pop) in the rest 
frame of the Higgs boson, which is shown in Fig. [T3l(a). The two Z bosons are reconstructed 
as explained above. Since the angle 6*^ and (pnp are defined in the rest frame of the Higgs 
boson, the non-zero transverse momentum of the Higgs boson does not affect these two 
variables. Therefore, as clearly shown in the figures, all the distributions of the angular 
variables mentioned above are the same for the RES, NLO and LO calculations. 
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Figure 13: Normalized distributions of coscpDP and cosOzz in the rest frame of the Higgs boson 
with mass 200 GeV in gg ^ H ^ ZZ ^ at the LHC. 



B. gg ^ H ^ ZZ ^ £+£~ui> 

Although the "gold-plated" mode, H ZZ , is considered to be the most 

effective channel for the SM Higgs boson discovery at the LHC, it suffers from the small 
decay branching of Z ^ . Moreover, the larger the Higgs mass becomes, the smaller 
the production rate is. When the Higgs boson mass is larger than 600 GeV, the H — > 
ZZ i^i'i^u channel may become important because the decay branching ratio (Br) of 
H ^ ZZ ^ is six times of the BioiH ^ ZZ ^ i+£~ The drawback is that 

one cannot reconstruct the Higgs mass from the final state particles due to the presence of 
two neutrinos. In this discovery channel, the missing transverse energy (^r) is crucial to 



suppress the background [37||. The distribution is shown in Fig. [T4l(a) which exhibits a 
Jacobian peak around M^/^, and the soft-gluon resummation effects smear the Jacobian 
peak and shift more events to the larger region. Similar to the H WW mode, the 
kinematics of this channel is similar to the W boson production and decay in the Drell-Yan 
process, therefore the shape of -^t distribution change significantly by the RES contributions. 
The Higgs boson mass can be measured from the peaks of the distributions of the transverse 
mass Mt and the cluster mass Mc, cf. Eq. ([3]), as shown in Fig. [T4l(b) and (c). Although 
the upper endpoint of is insensitive to high order corrections as we mentioned in the 
study of H ^ WW^*^ mode, the Jacobian peak is smeared out by the width {Th) effects 
of the Higgs boson. For Mh = 600 GeV, the total decay width of the Higgs boson is about 
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Figure 14: Normalized distributions of^T, and Mc in gg ^ H ^ ZZ £^£ uv channel with 
Mh = 600 GeV at the LHC. 

Table III: Acceptance of the process gg ^ H ^ ZZ ^ +f " for Mh = 140 (200) GeV and 

the process gg ^ H ^ ZZ — > i^i'vi? for Mh = 600 GeV after imposing cuts. 





Mh = 


140 GeV 


Mh = 


200 GeV 


Mh = 600 GeV 




basic (Eq. [8]) 


optimal (Eq. [9|) 


basic (Eq. E]) 


optimal (Eq. (H) 


basic (Eq. 


RES 


0.53 


0.15 


0.67 


0.14 


0.55 


NLO 


0.54 


0.12 


0.67 


0.11 


0.56 


LO 


0.53 





0.67 





0.58 



120 GeV, which is quite sizable and generates a noticeable smearing effect on the Jacobian 
peak. 

C. Acceptance study 



The discovery potential of the H ^ ZZ ^ £+£-£'+£'- and H ^ ZZ 



~vv modes 



has been studied in Ref. 



37l| after imposing the following cuts: 



For Mh = 140 GeV and 200 GeV, the intermediate mass range, we impose the basic 
cuts: 



> 7.0 GeV, \Yl\ < 2.5, > 20 GeV, 



and the optimal cuts: 



Pt > q ■ 



(8) 



(9) 
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where p§ and Yl are the transverse momentum and rapidity of each charged lepton, 
respectively, and p^^^"'' is the pr of the harder Z boson. 

• For Mh = 600 GeV, we require: 

p^>AOGeV, |Yl|<2.5, 
pif^ > 200 GeV, > 150 GeV, (10) 

where pff^ is the transverse momentum of the two charged lepton system. The numeri- 
cal results of the acceptances of the various cuts are summarized in Table lllli For the 
H ZZ^*^ mode, the RES and NLO contributions have almost the same 

acceptances after imposing the basic cuts. However, after imposing the optimal cuts the 
acceptance of the RES contribution is larger than the one of the NLO contribution by 25%, 
and the LO contribution is largely suppressed. For the H ZZ — > i^i^uu mode, the 
acceptances of the RES and NLO calculations are similar to each other. 



VI. CONCLUSION 



The search for the SM Higgs boson is one of the major goals of the high energy physics 
experiments at the LHC, and the vector boson decay modes, H WW^*^ or H ^ ZZ^*\ 
provide powerful and reliable discovery channels. The LHC has a great potential to dis- 
cover the Hi ggs boson even with low luminosity (~ 30fb~^) during the early years of run- 



nmg 
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38|, 



68|. In order to extract the signal from huge background events, we should 



have better theoretical predictions of the signal events as well as background events. In 
this paper, we examine the soft gluon resummation effects on the search of SM Higgs boson 
via the dominant production process gg ^ H at the LHC and discuss the impacts of the 
resummation effects on various kinematics variables which are relevant to the Higgs search. 
A comparison between the resummation effects and the NLO calculation is also presented. 

For H WW^*^ i+e-uu mode, we study Mh = 140 GeV at the Tevatron and 
Mh = 170 GeV at the LHC. Due to the spin correlations between the final state particles, 
this process is similar to the W boson production and decay in the Drell-Yan process. 
The shapes of the kinematics distributions are modified significantly by RES effects. For 
example, the effects could cause ~ 50% difference compared to NLO calculation in the 
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transverse momentum distribution of the leading lepton (p^^^'"'), when Mh = 170 GeV. The 
Higgs boson mass cannot be reconstructed directly from the final state particles because of 
two neutrinos. Therefore, the upper endpoint in the transverse mass distribution can be 
used to determine the mass of the Higgs boson, and we found that it is insensitive to the 
RES effects. After imposing various kinematics cuts, the LO, NLO and RES calculations 
yield similar acceptance of the signal events. 

For the H ZZ^*^ +f- mode, the so-called "gold-plated" mode, we study 

Mh = 140 GeV and Mh = 200 GeV at the LHC in this paper. We pursue an algorithm, 
called minimal deviation algorithm in this paper, to reconstruct the two Z bosons when the 
four charged leptons in the final state have the same fiavors. The RES effects change the 
shapes of kinematics significantly, e.g. p^"""" and distributions. However, the variables 
(f)DP and 9*^, defined in the Higgs rest frame, are insensitive to RES effects. After imposing 
the optimal kinematics cuts, the RES effects could increase the acceptance by 25% compared 
to that of NLO calculation while the LO contribution is largely suppressed. When the Higgs 
boson is heavy (600 GeV), we consider the H ZZ £^i~uu mode because of its larger 
decay branching ratio, as compared to the H ZZ — > i~^i~£'^£'^ mode. The shape of 
distribution, which is crucial to suppress the backgrounds, is largely modified because it is 
sensitive to the transverse momentum of the Higgs boson. 

In summary, we have presented a study of initial state soft-gluon resummation effects 
on the search for the SM Higgs boson via gluon-gluon fusion at the LHC. The effects not 
only significantly modify some of the kinematic distributions of the final state particles, as 
compared to the NLO and LO predictions, but also enhance the acceptance of the signal 
events after imposing the kinematic cuts to suppress the large background events. Therefore, 
we conclude that the initial state soft-gluon resummation effects should be taken into account 
as searching for the Higgs boson at the LHC. In addition, we note that the spin correlations 
among the final state leptons could be modified by the electroweak corrections to the Higgs 
boson decay. Therefore, we have implemented the NLO QED correction in the ResBos code, 
and the phenomenological study will be presented in the forthcoming paper. 
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